Finite set of invariants to characterize local Clifford equivalence of stabilizer states 
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The classification of stabilizer states under local Clifford (LC) equivalence is of particular impor- 
tance in quantum error-correction and measurement-based quantum computation. Two stabilizer 
states are called LC equivalent if there exists a local Clifford operation which maps the first state to 
the second. We present a finite set of invariants which completely characterizes the LC equivalence 
class of any stabilizer state. Our invariants have simple descriptions within the binary framework 
in which stabilizer states are usually described. 
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PACS numbers: 03.67.-a 

Stabilizer states constitute a class of multipartite pure 
quantum states which is of considerable interest in quan- 
tum information theory (QIT) and quantum computing 
(QC). A stabilizer state on n qubits is defined as a simul- 
taneous eigenvector of a maximal set of commuting ob- 
servables in the Pauli group, where the latter is the group 
generated by all n-fold tensor products of the Pauli matri- 
ces and the identity. To name but two of their presently 
most prominent applications, stabilizer states firstly ap- 
pear as codewords of quantum stabilizer codes in the the- 
ory of quantum error-correction and, secondly, they 
are the resources for QC in a computation model where 
only single-qubit measurements are performed (the one- 
way quantum computer 0). In order to understand the 
role of stabilizer states in these and other QIT and QC 
tasks, the mathematical and physical properties of these 
states have recently been studied in a number of accounts 
0, H H H 0, H One central problem in this research is 
the classification of stabilizer states under local unitary 
(LU) equivalence. Indeed, next to the natural relevance 
of this issue in the study of the entanglement properties of 
stabilizer states, it is also of particular importance both 
in the development of the one-way quantum computer 
and in the coding theoretic aspect of stabilizer states. In 
this paper we consider a restricted version of LU equiv- 
alence of stabilizer states, in that we consider only those 
local unitary operations which belong to the local Clifford 
group (LC). The local Clifford group consists of all local 
unitary operators which map the Pauli group to itself 
under conjugation. We will call two stabilizer states LC 
equivalent if there exists a local Clifford operator which 
maps the first state to the second. 

The restriction of considering LC equivalence rather 
than general LU equivalence has two main motivations. 
Firstly, the problem of general LU equivalence is difficult 
and it is therefore appropriate to first consider a more 
manageable subproblem, for which LC equivalence is a 
suitable candidate - given the very explicit connections 
between stabilizers states, the Pauli group and the (lo- 
cal) Clifford group. Secondly, this so-called 'subproblem' 
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might in fact not be a subproblem at all, as the natural 
question exists whether every two LU equivalent stabi- 
lizer states are necessarily LC equivalent. Although this 
question is to date unanswered in gene ral, our work j^], 
following work performed in Ref. |Toj , shows that this 
assertion is at least true for some interesting subclasses 
of stabilizer states and a positive answer to this question 
is widely believed. 

In the following we study LC invariants of stabilizer 
states, i.e. functions in the entries of a stabilizer state 
which take on equal values in LC equivalent states. We 
construct a finite complete set of invariants, which char- 
acterizes the LC equivalence class of any stabilizer state. 
Such a study of LC invariants primarily serves to obtain 
a characterization of LC equivalence classes of stabilizer 
states and to gain insight in the structure of these classes, 
and to a lesser extent to obtain algorithms which can 
recognize LC equivalence of two given states. Indeed, 
an efficient algorithm to recognize LC equivalence of sta- 
bilizer states has recently been presented in our earlier 
work and there is little hope that equally efficient algo- 
rithms can be constructed on the basis of invariants (see 
below). Therefore, we believe that the main merit of the 
present result is, firstly, that it gives a finite characteri- 
zation of the LC equivalence class of any stabilizer state 
and, secondly, that our invariants exhibit a very trans- 
parent structure in terms of the stabilizer formalism. 

Let us now start by recalling the basic notions con- 
cerning stabilizer states and the local Clifford group. A 
stabilizer state \tp) on n qubits is the unique simultaneous 
eigenvector with eigenvalue 1 of a set of n commuting and 
independent observables in the Pauli group Q n . The lat- 
ter consists of all operators of the form acr ai (g) • • • <g> o~ an , 
where a S {±l,±i} is an overall phase factor and a ai 
is either the 2x2 identity matrix or one of the Pauli 
matrices a x , o~ y , a z , for every i = 1, . . . ,n. The stabi- 
lizer iS(i/0 of is the set of all Pauli operators which 
have \tp) as an ei gen vector with eigenvalue 1. It is well 
known (see e.g. [ll|) that the stabilizer formalism has 
an equivalent formulation in terms of algebra over the 
field F2 = GF(2), where arithmetic is performed modulo 
2. The heart of this binary representation is an encoding 
of the Pauli matrices by pairs of bits, writing uq = 000, 
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a x — (Toij — (Tio and o y = on and subsequently 



(1) 



where u> := . . . , io n ), to' := (ty^, . . . , w' n ) G F£ 
are n-dimensional binary vectors. Note that the infor- 
mation about the ith qubit is distributed over the ith 
components of the vectors w and w' . Identifying every 
Pauli operator in an n-qubit stabilizer S(ip) with its 2n- 
dimensional binary index vector, one can show that S(tjj) 
corresponds to an n-dimensional self-dual linear subspace 
Sj, of F| n ; the self-duality of this subspace is with respect 
to a symplectic inner product on F|™. 

The Clifford group C\ on one qubit is the group of all 
2x2 unitary operators which map a u to a u a„r u -\ under 
conjugation, where u = x,y, z, for some a u — ±1 and 
some permutation 7r of {x, y, z}. The local Clifford group 
C l n on n qubits is the n-fold tensor product of C\ with 
itself. In the binary stabilizer framework, local Clifford 
operations U G C n correspond to nonsingular 2n x 2n 
binary matrices Q of the block form 



Q = 



A B 
C D 



(2) 



where the n x n blocks A,B,C,D are diagonal [jj. We 
denote the diagonal entries of A, B, C, D by <Zj, bi, Ci, di, 
respectively. The n submatrices 



Qi 



b, 

di 



S GL(2,¥ 2 ) 



(3) 



correspond to the tensor factors of U . We denote the 
group of all such Q by C l n . In the binary stabilizer frame- 
work, two n-qubit stabilizer states \<ip), \ip') are LC equiv- 
alent if and only if there exists an operator Q G C l n such 
that QSip — S^>, i.e., Q maps the space to the space 

Finally, let us recall the definition of support in this 
context of stabilizers, which will be of considerable im- 
portance in the following. The support supp(w) of any 
vector v G F| n is the set 

supp(u) = {i G {1, . . . ,n} | (vi,v n+i ) / (0,0)}. 

Note that supp(w) contains exactly those i G {1, . . . , n} 
such that the ith tensor factor of a v differs from the iden- 
tity. 

We are now in a position to state the first result of 
this paper. In theorem 1 a finite set of invariants which 
characterizes the LC equivalence class of any stabilizer 
state, is presented. 

Theorem 1 Let \ip) be a stabilizer state on n qubits. 



Let r G No and consider subsets ur 



C {!,..., n} 



for every k,l G {l,...,r} with k < I. Denote Q := 
(w 1 , uj 2 , . . . , lu 12 , ui 13 , . . . ) and let T^ r {tp) be the set con- 
sisting of all tuples (v , ■ ■ . , v r ) G S^, x • • • x S 1 ^, satisfying 



SUpp(v k ) = Ld k , SUpp(v k 



(4) 



Then (i) |7^p r (-0)| is an LC invariant and (ii) the LC 
equivalence class of \tp) is completely determined by the 
values of all invariants j7^p ra (^)| (i.e. where r = n). 

Proof: Statement (i) is trivial. To prove (ii), let \ip) 
and be two stabilizer states on n qubits such that 
l^nWl = l^nW)! f° r every tuple Q as defined above. 
We will show that this implies that \ip) and are LC 
equivalent. Fix a basis {v 1 , . . . , v n } of and define the 
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tuple of sets = (u>l, . . . , u>. 

uj k = supp(w fc ), Lu kl = supp(v k 



such that 



(5) 



for every k,l — l,...,n with k < I. It follows that 
l^rPnWI i s nonzero, and therefore \T^{\p')\ is nonzero 
as well, since these numbers are equal by assumption. 
Therefore, there exist vectors w 1 , . . . ,w n G S^i such that 



supp(u; fc ) = uj k , supp(w k +w l )= Lu kl , 
for every k,l. Expressions (J3J) and JHJl show that 

K fe ,^ + . i ) = (o,o) iff K fc ,^ +4 ) = (o,o) 

(v^Vn+i) = U ,V l n+i ) iff ( 



(6) 



for every i, k, I — 1, ... ,n. This implies the existence of 
n operators Qi G GL(2,¥2) such that 



Qi 



v k 




w k 


_ V n+i _ 







(7) 



for every i, k = l,...,n, and the operators Qi con- 
sequently constitute an operator Q G C l n such that 



Qv = w for every k = 1, 



Therefore, the space 



QSip is spanned by the vectors w . . ,w n . Moreover, 
these vectors are linearly independent, as they are the 
images of a linearly independent set under an invertible 
linear transformation. Therefore, {w , . . . , w n } is a basis 
of S^i (since the latter is an n-dimensional vector space) 
and we obtain QS^p = S^p> . This proves the theorem. □ 
Let us briefly discuss this result. The invariants in 
theorem 1 are presented as cardinalities of certain subsets 
of $0 x • • • x , which are defined in terms of simple 
constraints on the supports of their elements. Note that 
for r — 1 these invariants simply count the number of 
elements in the stabilizer with a prescribed support, i.e., 
for every lo C {1, . . . , n} one has an invariant 



|{« £ S^i | SUpp(f) = L0}\ 



(8) 



These invariants are in fact 'local versions' of the so-called 
weight distribution of a stabilizer, a well known concept 
in quantum (and classical) coding theory which is e.g. 
used to produce bounds on how good codes can be (see 
e.g. [12 )■ Our invariants for r > 2 are in this respect 
generalizations of this well known notion. 

While the structure of the invariants in theorem 1 is 
indeed very transparent, a direct calculation of the num- 
bers \T^ r (ip)\ is likely to be hard. It is in this context 
interesting to consider the following variant of theorem 
1, which presents an alternative complete family of LC 
invariants: 
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Theorem 2 Using the same notations as in theorem 1, 
let V„ r (^>) be the vector space consisting of all tuples 
(v v r ) £ St£ x • ■ • x satisfying 



supp(v k ) C uj k , supp(v k + v l ) C LO k 



(9) 



T/ien (zj t/ie dimension of V^ r {if>) is an LC invariant 
and (ii) the LC equivalence class of\i/)} is completely de- 
termined by the values of all invariants dim(V^ n (V0) • 

Proof: Statement (i) is again trivial. We now prove 
(ii). Let lu C {1, . . . , n} and define the functions S u , e w : 
F 2 " — ► C, respectively, by 6u(w) = 1 if supp(w;) = to and 
<5w(w) = otherwise, and e u (w) = 1 if supp(w) C u> and 
£uj{w) = otherwise. It is then straightforward to show 
the following relations: 



U)' C-LJ 

s* = (-i) H E ^ 



(10) 



the first of which is trivial and the second of which can be 
verified using straightforward combinatorics. Letting f2 
be a tuple of sets uj k , co kl C {1, . . . , n}, with k, I = 1, . . . , r 
and k < I, it follows that 




iv" (VOI 



where the sums run over all r-tuples (v , . . . , v r ) £ S^ r . 
Using (fTU|) one finds that every invariant |^p°(-)| is a lin- 
ear combination of the invariants {|V^ 2 r .(-)|}a with con- 
stant coefficients (and vice versa). The proof then read- 
ily follows after applying theorem 1 and noting that the 
dimension of any vector space V over F2 is given by 
dim(V) =log 2 |V|. □ 
The main advantage of the invariants dim(V^ 2 r ( - 0)) 
over the \%^ r (ip)\'s is that in the former case one is 
dealing with (dimensions of) vector spaces, which are 
much more manageable from a computational point of 
view. Let us illustrate this with a simple example: take 
u> C {1, . . . , n} as before and consider the invariant 



\{v £ I supp(w) C uj}\ 



(11) 



which is the counterpart of jHJ in theorem 2. This in- 
variant can indeed be calculated efficiently as follows: let 
v 1 , . . . ,v n be a basis of S^, and let S = [v 1 ] . . . \v n ] be 
the In x n binary matrix which has the vectors v J as its 
columns (the matrix S is generally referred to as a gener- 
ator matrix of the stabilizer). Denoting by Sf the 2 x n 
submatrix of 5* which is obtained by assembling the ith 
and the (n+i)th row of S 1 , for every i — 1, . . . , n, one read- 
ily verifies that the invariant (|llfl is equal to the corank 



(i.e. the dimension of the kernel) over F2 of the matrix 
Su> '■= where Q denotes the complement of a; in 

{1, . . . , n}. Note that the calculation of the (co)rank of a 
matrix over any field indeed requires only a polynomial 
number of operations in the dimensions of this matrix. 
This shows that the invariant ljll|l can easily be calcu- 
lated when a generator matrix of the stabilizer is known 
(note that a stabilizer state in indeed usually presented in 
terms of a generator matrix). A straightforward contin- 
uation of the above argument shows that, e.g. for r = 2, 
an invariant dim(V£ 2 2 ('i/>)); where = (w 1 , uj 2 , uj 12 ) with 
uj 1 , lu 2 , lu 12 C {l,...,n}, is equal to the corank of the 
matrix 



(12) 



and analogous results hold for arbitrary degrees. This 
shows that the invariants dim(V^ r (^>)) can indeed be 
evaluated efficiently; moreover, the information needed 
to perform the calculation (i.e., the matrices 5 U ) can eas- 
ily be extracted from a generator matrix of the stabilizer. 

There is in fact a close relationship between the LC in- 
variants dim(V^ r (V>)) and the polynomial LU invariants 
that we considered in 6]. In both cases, invariants corre- 
spond to dimensions of certain subspaces of X • • • X S^, 
that are defined in terms of constraints on the supports 
of their elements. In fact, if r equals 1 or 2 then the 
LC invariants in theorems 1 and 2 are also LU invari- 
ants; in other words, the invariants dim(V^ r (V>)) and 
7^,(^)1 take on equal values on LU equivalent stabi- 
lizer states. This property is well known (see e.g. [la ]') 
for the invariants JSJ and i|ll|) . corresponding to r = 1. 
In order to prove this assertion for r = 2, let W = 
be (the projector associated with) a stabilizer state and 
take il = (uj 1 , uj 2 , uj 12 ) . Then one can readily verify that 
l^n 2(^)1 iS: U P t° a multiplicative constant, equal to 



Tr {(Tr^M/) (Tv^) (Tr^)} , 



(13) 



which is manifestly invariant under the action of local 
unitary operators (the operation Tr w denotes tracing out 
all qubits in the set uj). Moreover, this also shows that 
every I^^WOI is an LU invariant, as it can be writ- 
ten as a linear combination with constant coefficients 
of the |V^ 2 (V0I' S - Summarizing the current discussion, 
if r equals 1 or 2 then the invariants dim(V„ r (ip)) and 
\Tnr (VOI t a k e 011 e< l ua l values on LU equivalent stabilizer 
states. It is clear that this result is very interesting in 
relation with the problem whether LC equivalence and 
LU equivalence of stabilizer states are identical notions. 

Finally, we wish to note that the two main results in 
this paper are to be regarded as upper bounds on the 
maximal r which needs to be considered in order to ob- 
tain a complete set of invariants. Indeed, it is likely that 
only invariants corresponding to (substantially) lower r 
need to be considered in order to recognize LC equiva- 
lence between two stabilizer states. One indication for 
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this conjecture is the proof of theorem 1, which shows 
that all the information of the LC equivalence class of a 
given state \ip) is in fact contained within a single invari- 
ant |7^p*(^)|. A second state \ip') is LC equivalent to \ip) 
if and only if 

Kn(i>)\ = \X?;M% (i4) 

(in fact, even \T^(yj')\ > is sufficient) and the infor- 
mation in all other invariants is in fact redundant. It 
is not unlikely that there exist smaller complete lists of 
invariants (i.e. of smaller r) which exhibit less redundan- 
cies. Secondly, our work shows examples of LC equiva- 
lence classes which are indeed characterized by invariants 
of small r: e.g., any stabilizer state is LC equivalent to 
the GHZ state on n qubits if and only if the values of 
the invariants coincide for both states 0- We are 
therefore led to believe that the results in theorems 1 
and 2 can in principle be improved - if not for all stabi- 
lizer states then at least for some interesting subclasses of 
states. It is in this context interesting to point out that 
it has been conjectured and later disproved in a graph 
theoretical context that the list of LC invariants are 
sufficient to characterize the LC equivalence class of all 
stabilizer states |5. Il l|. Therefore, invariants where r is 
at least 2 must be taken into account. It is to date not 
clear whether also higher r are to be considered. 

Summarizing this paper, we have characterized LC 
equivalence classes of stabilizer states by means of a finite 



set of invariants. Our invariants have simple descriptions 
within the binary stabilizer framework, in that they are 
equal to dimensions of certain subspaces of x ■ • • x S^. 
Furthermore, we discussed the link with existing families 
of LU invariants. 
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